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Let G be a graph on n vertices with r := [n/2j and let Ai > • • • > A„ be adjacency 
eigenvalues of G. Then the Hiickel energy of G, HE(G), is defined as 



The concept of Hiickel energy was introduced by Coulson as it gives a good approximation for 
the 7r-electron energy of molecular graphs. We obtain two upper bounds and a lower bound 
for HE(G). When n is even, it is shown that equality holds in both upper bounds if and only 
if G is a strongly regular graph with parameters (n, k, A, /i) = (4f^ + 4t + 2, 2t'^ + 3t + 1, + 
2t, + 2t + 1), for positive integer t. Furthermore, we will give an infinite family of these 
strongly regular graph whose construction was communicated by Willem Haemers to us. He 
attributes the construction to J.J. Seidel. 
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1 Introduction 

Throughout this paper, all graphs are simple and undirected. Let G be a graph with n vertices 
and A be the adjacency matrix of G. Then the eigenvalues of G are defined as the eigenvalues of 
A. As all eigenvalues of A are real, we can rearrange them as Ai > • • • > A„. I. Gutman (see [9]) 
defined the energy of G, E(G), by 



In chemistry, the energy of a given molecular graph is of interest since it can be related to the total 
TT-electron energy of the molecule represented by that graph. The reason for Gutman's definition 
is that E(G) gives a good approximation for the 7r-electron energy of a molecule where G is then 
the corresponding molecular graph. For a survey on the energy of graphs, see |9]. The Hiickel 
energy of G, denoted by HE(G), is defined as 



The idea of introducing Hiickel energy (implicitly) exists in Erich Hiickel's first paper [10 in 1931 
and is also found in his book [11]. The concept was explicitly used in 1940 by Charles Coulson [2] 
but, most probably, can be found also in his earlier articles. In a "canonical" form, the theory 
behind the Hiickel energy was formulated in a series of papers by Coulson and Longuet-Higgins, 
of which the first (and most important) is [3] . In comparison with energy, the Hiickel energy of a 
graph gives a better approximation for the total 7r-electron energy of the molecule represented by 
that graph, see |7]. Clearly for a graph on n vertices, HE(G) < E(G), and if G is bipartite, then 
equality holds. Koolen and Moulton in |12l I13j gave upper bounds on the energy of graphs and 
bipartite graphs, respectively. These bounds have been generalized in several ways. Obviously, the 
upper bounds of Koolen and Moulton also give upper bounds for the Hiickel energy of graphs. In 
this paper, we obtain better upper bounds for the Hiickel energy of a graph. More precisely, we 
prove that for a graph G with n vertices and m edges, 



n 



E(G) :=5^|A,|. 





(1) 



if n is even, and 





if n is odd. Then we show that 



HE(G) < 2 (1 + V^^) , 



(3) 
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if n is even, and 

RE{G)<^(^i+v^-^y (4) 

if n is odd. Moreover, equality is attained in ([l]) if and only if equality attained in ([s]) if and only if G 
is astrongly regular graph with parameters (n, k, A, ^) — (4t^+4t+2, 2t^+3i+l, t'^+2t, <^+2i+l). 
The proofs of the above upper bounds are given in Section 2. It is known that E(G) > — 1 
for any graph G on n vertices with no isolated vertices with equality if and only if G is the star 
Ki,n-i (see [4 ). In Section 3, we prove that the same bound holds for Hiickel energy. In the last 
section, we give a construction of srg(4i2 +4t + 2, 2t^ + 3t + 1, + 2t, f + 2t + 1). 



2 The upper bound for Hiickel energy 

In this section we prove ([TJ, and Q. The equality cases are also discussed. We begin by 

stating a lemma which will be used later. 



Lemma 1. Let G be a graph with n vertices and m edges. Suppose r ln/2\ and 

r 



i=l 

If m > n — I > 2, then 



a Aw} 

- < — • 5 



Proof. First, suppose m > n. Then G contains a cycle, and so by interlacing, we see ^n + ^n-i > 
2. Therefore, a/r < {2m— 2) /r < Arri^ /n^ . If m = n — 1 and G is connected, then G is a tree. Thus 
a — m, and obviously (|5| holds. So in the rest of proof we assume that G is disconnected and 
m — n~\. If G has at least three two non-trivial components, then at least one of the components 
contains a cycle. The component containing a cycle has either an eigenvalue at most < —2, or 
two eigenvalues < — 1 and the other two components have an eigenvalue < — 1. It turns out that 
+ ^n-i + + > 4 whcrc n > 7. Hcncc 

a < 2m - 4. (6) 

It is easily seen that ([5| follows from ([6|. Now, suppose that G has two non-trivial connected 
components Gi and G2, say. Let Gi,G2 have ni,n2 vertices and mi,m2 edges, respectively. 
First suppose ni,n2 > 3. If Gi or G2 contains a Ki^2 as an induced subgraph, we are done by 
interlacing. So one may assume that both Gi and G2 contain a triangle. It turns out that mi > ni 
and TO2 > n.2. Hence G must have an isolated vertex which implies n > 7. On the other hand, 
by interlacing, the four smallest eigenvalues of G are at most —1 implying M. Now, assume that 
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Figure 1: The paw and the diamond graphs 



ni — 2. So G2 must contain a cycle Ci. We may assume that Cg has no chord. If ^ > 4, we are 
done. So let £ = 3. If 712 = 3, then G does not have isolated vertices, i.e., G — C3U K2, for which 
([5]) holds. Thus ri2 > 4 which means that at least one of the diamond graph, the paw graph (see 
Figure [ij, or K/^ is induced subgraph of G. If it contains either the diamond or the paw graph, 
we are done by interlacing. If it contains K4, then G must have at least two isolated vertices, i.e., 
n > 8. Thus the four smallest eigenvalues of G are at most —1 which implies Finally, assume 
that G has exactly one non-trivial component Gi with ni vertices. It turns out that ni > 3 and 
Gi contains a cycle. By looking at the table of graph spectra of [5l pp. 272-3], it is seen that if 
Til =3,4, G satisfies ([5]). If ni > 5, then making use of the table mentioned above and interlacing 
it follows that Xf^ + X^^i + — 4 unless Gi is a complete graph in which case G has at least 
11 vertices and the four smallest eigenvalues of G is —1, implying the result. □ 

2.1 Even order graphs 

In this subsection we prove ([ij and ([3| for graphs of an even order. The cases of equalities are also 
characterized. 

Theorem 2. Let G be a graph on n vertices and m edges where n is even. Then |7p holds. 
Moreover, equality holds if and only if n = 4t^ + 4t + 2 for some positive integer t, m = {2t^ + 2^-1- 
l){2t^ + 3< + 1) and G is a strongly regular graph with parameters (n, k, A, ji) = (4i^ + 4i + 2, 2t^ + 
3t+l, t^ + 2t, <2 + 2t+l). 

Proof. Let n = 2r and Ai > A2 > • • • > A„ be the eigenvalues of G. Then 



^ A, = and ^ \l = 2m. 




Let a be as in Lemma [T] Then 



n 
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and Af < a < 2m — 2. By the Cauchy-Schwartz inequality, 



HE(G) = 2^A, < 2Ai +2y'(r- l)(a- A?). 
1=1 

The function x i-^ x + \J {r — 1)(q; — x^) decreases on the interval \J cajr < x < ^/a. By Lemma [l| 
m/r > y^a/r. Since Ai >m/r, we see that 

HE(G) < /i(a) — + 2^{r ~ l){a - /r^). (7) 

On the other hand, 

n 

HE(G) = -2 ^ A, < /2(a) := 2^r{2m-a). (8) 

Let 

f{a) := min{/i(Q:),/2(a)}. 

We determine the maximum of /. We observe that /i and /2 are increasing and decreasing functions 
in a, respectively. Therefore, max / — f{ao) where oq is the unique point with /i(ao) = /2('^o)- 
So we find the solution of the equation fi{a) — /2(a)- To do so, let a :— \/ a ~ /r^ and consider 
the equation 

m 
r 

This equation has the roots 



h (T\/r — 1 = -^Z r{2m — cr^ — m^/r^). 

-m\/r — 1 ± ry^2'm(2r^ — r — m) 



o'i,2 : 



r(2r- 1) 

If m < 2r^/{r + 1), then (Ti > and so 

2m 2^r - 1 



HE(G) < — + — ^-^ f-mVr - 1 + rJ2m(2r'^ - r - m)) 

r r[2r — 1) V / 

2m y^2m{n — 2) (n^ — n — 2m) 



n — 1 n — 1 

otherwise 



(9) 



HE(G) < 2\/r(2m-m2/r2), (10) 

which is less than 4m /n for m > 2r^/{r + 1). This shows that Inequality ([l]) holds. 

Now let us consider the case that equality is attained in (jlj. First let m < ^q^. Then equality 
holds if and only if 

2. X2 = ■ ■ ■ — Xr = 77=1' 
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3. \r+i = • • • = A„ = -^^2™ - CTi - ni^/r"^. 

The first condition shows that G must be "-regular, and the second and third conditions imply 
that G must be strongly regular graph as a regular graph with at most three distinct eigenvalues 
is strongly regular, cf. [8, Lemma 10.2.1]. From (S) Lemma 10.3.5] it follows that G has to have 
the parameters as required in the theorem. If m > ^q^, then with the same reasoning as above 
one can show that G has to be strongly regular graph with eigenvalue of multiplicity r — \ and 
by [HI Lemma 10.3.5] such a graph does not exist. □ 

Optimizing the Hiickel energy over the number of edges we obtain: 

Theorem 3. Let G he a graph on n vertices where n is even. Then ^ holds. Equality holds if and 
only if G is a strongly regular graph with parameters (4t^ + At + 2, 2t^ + 3i + 1, + 2t, + 2t + 1), 
for some positive integer t. 

Proof. Suppose that G is a graph with n vertices and m edges. If m < n — 2, then ^ obviously 
holds as E(G) < 2m (see |9]). If m > n — 1, then using routine calculus, it is seen that the right 
hand side of Q — considered as a function of m — is maximized when 

m — n{n — 1 + \/ n — 1) /4. 

Inequality ([3| now follows by substituting this value of m into ([ij. (We note that the maximum 



of the right hand side of (10) is 2n^/{n + 2) which is less than the above maximum.) Moreover, 
from Theorem [2] it follows that equality holds in (|3| if and only if G is a strongly regular graph 
with parameters (4*^ + U + 2, 21"^ + 3t + I, + 2t, t'^ + 2t + I). □ 

2.2 Odd order graphs 

In this subsection we prove ([ij and ^ for graphs of an odd order and discuss the equality case 
and tightness of the bounds. 

Theorem 4. Let m > n — 1 > 3 and G be a graph with n vertices and m edges where n is odd. 
Then holds. 

Proof. Let n = 2r + 1, a be as before, and 
We have 

2m- a > (r + l)/?^. 
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This obviously holds if /? < 0. For /3 > it follows from the following: 

2ni-a-B'= > ' A? > - I > ' A,-^ 



i—r-\-2 \i—r-\- 



=r+2 

2 

^ 2 a2 



where the first inequality follows from the Cauchy-Schwartz inequality. In a similar manner as the 
proof of Theorem[2] we find that HE(G) < min{/i(a, /3), /2(q;,/3)}, where 



13) = 4m/n + 2v/(r' - 1) (a - 4m2/n2) + /?, and (11) 
/2(a,/3) = 2v/r(2m-a-/32)-/?. (12) 

Let 

/(a,/3) :=min{/i(a,/3),/2(a,/3)}. 
We determine the maximum of / over the compact set 

D := {{a,f3) : a> Am'^/n'^, 2m- (r + l)/3^ > a} . 

Note that for (a, (3) ^ D one has — /3o < P < Po, where 



2 /m(n2 - 2m) 
n \ n + 1 

Neither the gradient of /i nor that of /2 has a zero in interior of D. So the maximum of / occurs 
in the set 

L:={(a,/3):/i(a,/3) = /2(«,/3)}, 
where the gradient of / does not exist or it occurs in the boundary of D consisting of 

Di := {(a,/3) : a = W/n', ~Po < l3 < /3o}, 

D2 := {(a, p):a = 2m-ir + l)/?^, -/3o < /3 < /3o}. 

For any (a, /3) G _D, /2(a, /3) < f2{^m'^/n^i0)- It is easily seen that the maximum of f 2(4:771^ /n^ , /3) 
occurs in 

1 / 2m{n? — 2m) 



Pi — u 

n V 2n — 1 
Therefore, 

max /2 = /2(4TOVr^^/3l) = -v/2m(2n - l)(n2 - 2m). (13) 
In the rest of proof, we determine max / for 

m < r; (14) 

- 2(n2 -4n+ 11 ' ^ ^ 
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if ™ > 2(k^^4«+ii) ' '2 1 follows from (jisj). 
On Z?i, wc have 



max /|^^ < /i(4mVn2, /3o) = 4m/n + /Jq. 



On Z?Q, one has 



= 4m/n + 2^(r - l)(2m - (r + l)/?^ - 4m2/n2) + /?, and 
/2(/3)-(n-l)|/3|-/3. 

In order to find max , we look for the points where = /2(/3). The solution of this equation 

for ^ < is 

^ -~2m{n + 1) — yj2Tan(n^ — 2n — 3)(n2 — n — 2m) 

n(n2 — 1) ' 

and for /3 > is 

2m(n — 3) + y^2m(n — 3){n'^ — An^ — 2mii? + 'drfi + 6mn — 8m) 

= zrr^ 



n{'n? — 4n + 3) 

We have (32 > -/3o if and only if m < "2^^n+3} ' ^^'^ < ;9o if and only if m < 2(n2'-4n+ii) ■ 
Moreover ^(/Ja) > ^(/Sa). Thus if m < 2(^^^^^i^, 

, , 2m(n + 1) + j2mn{Ti? - 2n - 3)(n2 - n - 2m) 
max = /2(/32) = ^ ^^3^ . 



Now we examine max/|^. Let a := \J a — Am? /v? . To determine (a,/3) satisfying fi{a,P) = 
/2(a, /3) it is enough to find the zeros of the following quadratic form: 

(2n - 4)cr2 + 4(2m/n + /3)V^r^cr + (4m/n + 2(3f -{n- l)(4m - 2/3^ - 8m^/n^). (15) 

The zeros are 



— (2m + n[3)^/n — 3 ± ■y/(n — l)(2mn3 — + /J^n^ — Amn? — Amfin ~ Anm? + 4m2) 

n(n — z) 



Note that CT2 < and so is not feasible. Let us denote the constant term of (15 1 by h{(3) as a 
function of -/3o < fi < (3a- Then cti > if and only if h{(i) < 0. Moreover /i(/3) < h{l3o), and 
hWo) < if 

n2(7i-3)2 
- 2(n2-4n+ll)' 

Thus, with this condition on m, /i becomes /i(/3) = 4m/n + 2y/^F^ai. The roots of /((/3) = 
are 

-2m(n2 - 3n + 1) ± ^2mn(n2 - 3n + l)(n2 _ ,7 _ 2m) 

P4,5 = 



i{n^ - An? + 4n - 1) 
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It is seen that —(3o < /^s < /?4 < unless m < n^/(2n — 2). We have 

2y/2mn{n'^ — 3n + l){n'^ — n — 2m) 



i{n - 2)(7i3 - An? + 4n - 1) 



Thus fi^Pi) > fiiPs)- It turns out that fi decreases for P > P4, so /i(/?4) > /i(/3o)- It is easily 

2(n2_4„+ii) 



seen that /i(/34) > /i(-/5o)- Therefore, for m < 2(n^'^-4n+ii) ^^'^^ 



./ON 2m ^2mn(n2 - 3n + l){'n? - n - 2m) 
max /|^ = /i(/34 = T + TT . (16) 

The result follows from comparing the three maxima max /|^^ , max /|^^ , and max /|^ . □ 
Theorem 5. Let G be a graph on n vertices where n is odd. Then holds. 



(17) 



Proof. The maximum of the right hand side of ( |13[ ) — as a function of m — is 

n(n — 3)V2(n+l)(2n- 1) 
n2-4n + ll 



which is obtained when m is given by (14 1. Also, the right hand side of (16 1 is maximized when 

n(n — 1 + ^/n) 



This maximum value is equal to ^ \^ ~^ ^ ^ j which is greater than ( 17 1. This completes the 
proof. 



□ 

Remark 6. Here we show that no graph can attain the bound in ([2|. Let us keep the notation 
of the proof of Theorem|4] First we consider m > ^n^{n — 3)^/ (n^ ~ An+ 11). Therefore, HE(G) 
equals (13). Then a = im'^/n'^ and A^+i = (3i. This means that G is a regular graph with only one 
positive eigenvalue. Then by [5] Theorem 6.7], G is a complete multipartite graph. As the rank 
of a complete multipartite graph equals the number of its parts, G must have r + 2 parts. Such a 
graph cannot be regular, a contradiction. Now, we consider to < ln^(n — 3)^/(?i^ — 471 + 11). Hence 



HE(G) is equal to (16 1. Then G must be a regular graph of degree k, say, with A2 = • • • = A^, 
Ar+i — Pa, and \r+2 = ■ ■ ■ — A„. Since A^+i — Pa < 0, A2 and A„ have different multiplicities, and 
thus all eigenvalues of G are integral. Let A2 = t. Then A„ = — t — s, for some integer s > 0. It 
follows that k + A^+i = t + rs. This implies that either s = or s = 1. If s = 0, then k + A^+i = t. 



and so HE(G) = fc + (n — 2)t. This must be equal to ( 16 1 which implies 



k + ^y nk{n? — 3n + l)(n — 1 — fc) 
^ n2 - 3n + 2 ' 

Substituting this value of t in the equation nk = k^ + {n — 2)t^ + {t — k)"^ and solving in terms 
of k yields fc = n/(n — 1) which is impossible. If s = 1, then k + Xr+i = t + r, and so HE(G) = 
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k + {n-2)t + r. It follows that 

_k-{n- 1)2/2 + ^nk(n^ - Sn + l)(n - 1 - fc) 

Substituting this value of t in the equation nk ~ fc^ + (r — 1)^^ + r{t + 1)^ + (r + t — fc)^ and 
solving in terms of k yields k = (n — I + y'n) /2 which implies t = {\/n — l)/2. Therefore, we have 
A^+i = — a contradiction. 

Remark 7. Note that a conference strongly regular graph G, i.e, a srg(4i + 1, 2t, t — 1, t), 
has spectrum {[2t]\[{-l + y/UTl)/2]'^* ,[{-1 - V^tTl)/2]'^*} , and hence HE(G) = ^(1 + 
V4F+T). This is about half of the upper bound given in Q. For odd order graphs, we can come 
much closer to Q. Let G be a srg(4i2 +At + 2, 2t^ + 3t + 1, + 2t, t^ + 2t+\). If one adds a new 
vertex to G and join it to neighbors of some fixed vertex of G, then the resulting graph H has the 
spectrum 

where Ai, A2, A3 are the roots of the polynomial 

p{x) - {2f + ■it)x'^ - [bt^ + 7t + 2)x + + IQt^ + 8f + 2t. 

It turns out that p(-\/2<) = (4 + 'i^/2)t^ + (8 + ^^/2)t^ + (2 + 2^/2)t. Hence A3 < -^/2t and so 

HE(i?) > 2(2^2 + 2t)(t+l)+2y2t= '^^'^^ ^ ( ^4^^ + 4^ + 3) + 0(4^^ + 4t + 3). 
This shows that Q is asymptotically tight. 

3 Lower bound 

It is known that E(G) > 2\/n — 1 for any graph G on n vertices with no isolated vertices with 
equality if and only if G is the star i^i^„_i (see [1]). Below we show that this is also the case for 
Hiickel energy. 

Theorem 8. For any graph G on n vertices with no isolated vertices, 

HE(G) > 2\Ar^. 
Equality holds if and only if G is the star Ki ^-i- 

Proof. If Gi, G2 are two graphs with ni,n2 vertices, then HE(GiUG2) > HE(Gi) +HE(G2), and 
^Jni — 1 + ^/ ni — \ > -y/ni + 712 — 1 for ni, 712 > 2. This alows us to assume that G is connected. 
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The theorem clearly holds for n < 3, so suppose n > 4. Let p^q be the number of positive and 
negative eigenvalues of G, respectively. Let n be even; the theorem follows similarly for odd n. If 
p— I, then by [5J Theorem 6.7] , G is a complete multipartite graphs with s parts, say. If s < f + 1, 
HE(G) = E{G) and we are done, so let s > ^ + 2. Note that the complete graph Ks is an induced 
subgraph of G, so by interlacing, A„_s+2 < ^ L Therefore X^i < —1 and thus IIE(G) > n and this 
is greater than 2\/n ~ 1 for n > 3. So we may assume that p> 2. This implies that q >2 as well. 
Note that either p < ^ + 2 or q < f + 2. Suppose q < f + 2; the proof is similar for the other 
case. If we also have p < f + 2, we are done. So let p > f+2. We observe that 

E(G) = 2(Ai + • • • + Ap) 

< 2(Ai H h A^) + 2(Ap_2r H \- K) 



It turns out that 



< HE(G) + - — ^HE(G). 



HE(G) > I^E(G). 



On the other hand, we see that E(G) > p + 2, as the energy of any graph is at least the rank of its 
adjacency matrix ([6], see also [1]). Combining the above inequalities we find 

HE(G)>^(p + 2) 
n n 

> 



2(n - 2) ' 



and the last value is greater than 2\/n — 1 for n > 4. □ 



4 A construction of srg(4t2 + 4t+2, 2t^ + 3t + l, t^ + 2t, t^ + 2t+l) 

In this section we give an infinite family of strongly regular graphs with parameters {n,k,X, ij) = 
(4i^ + 4t + 2, 2i^ + 3t + 1, + 2t, + 2t + 1) whose construction was communicated by Willem 
Haemers to us. He attributes the construction to J.J. Seidel. 

Let G be a graph with vertex set X, and Y Q X. From G we obtain a new graph by leaving 
adjacency and non-adjacency inside Y and X \ y as it was, and interchanging adjacency and non- 
adjacency between Y and X\Y . This new graph is said to be obtained by Seidel switching with 
respect to the set of vertices Y . 

Let q = 2i + 1 be a prime power. Let F be the Paley graph of order q^ , that is, the graph 
with vertex set GF(g^) and x ^ y ii x — y is not a square in GF(g^). Let a; be a primitive element 
of GF(g2) and consider V = | i = 0, ...,(?- 1} U {0}. Then V is the subfield GF{q) 
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of GF(g^) and V forms a coclique of size q. Now {a + V \ a G GF{q'^)} forms a partition into 
q cocliques of size q. Add an isolated vertex and apply the Seidel switching with respect to the 
union of t disjoint cocliques. The resulting graph is a strongly regular graph with parameters 
(4t^ + 4t + 2, 2t^ + t, — 1, P). Taking the complement of this graph give us a strongly regular 
graph with the required parameters. Note that for t = I, there exists only one such a graph namely 
the Johnson graph J(5,2), for t = 2 there are 10. For t < 9 they do exist and t ^ 10 seems the 
smallest open case. 

Acknowledgements. We are very grateful for the discussions with Patrick Fowler and Ivan 
Gutman. Partick Fowler suggested the name of Hiickel energy and provided us with the reference 
[7] . Ivan Gutman presented us the history of the Hiickel energy and provided us with the references 
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